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ABSTRACT 

Observations suggest that the properties of Type la supernovae (SNe la) may depend on environmental 
characteristics, such as morphology, metallicity, and age of host galaxies. The influence of these envi- 
ronmental properties on the resulting SNe la is studied in this paper. First it is shown that the carbon 
mass fraction X{G) in the C-l-O white dwarf SN la progenitors tends to be smaller for lower metallicity 
and older the binary system age. It is then suggested that the variation of A(C) causes the diversity in 
the brightness of SNe la: a smaller X{C) leads to a dimmer SN la. Further studies of the propagation 
of the turbulent flame are necessary to confirm this relation. Our model for the SN la progenitors then 
predicts that when the progenitors belong to an older population or to a low metallicity environment, 
the number of bright SNe la is reduced, so that the variation in brightness among the SNe la is also 
smaller. Thus our model can explain why the mean SN la brightness and its dispersion depend on the 
morphology of the host galaxies and on the distance of the SN from the center of the galaxy. It is further 
predicted that at higher redshift (z > 1) both the the mean brightness of SNe la and its variation should 
be smaller in spiral galaxies than in elliptical galaxies. These variations are within the range observed in 
nearby SNe la. In so far as the variation in X{G) is the most important cause for the diversity among 
SNe la, the light curve shape method currently used to determine the absolute magnitude of SNe la can 
be applied also to high redshift SNe la. 

Subject headings: stars: white dwarfs — stars: evolution — binaries: close — supernovae: general — 
cosmology: miscellaneous 



1. INTRODUCTION 

It is widely accepted that Type la supernovae (SNe 
la) are thermonuclear explosions of accreting C+0 white 
dwarfs (WDs), although the nature of the progenitor bi- 
nary system and the detail of the explosion mechanism are 
still under debate. 

SNe la are good distance indicators, and provide a 
promising tool for determining cosmological parameters 
(e.g., Branch and Tammann 1992). From the observations 
of high redshift SNe la, both the SN Cosmology Project 
(Perlmutter et al. 1999) and the High-z SN Search Team 
(Riess et al. 1998) have suggested a statistically signifi- 
cant value for the cosmological constant. However, SNe la 
are not perfect standard candles, but show some intrinsic 
variations in brightness. When determining the absolute 
peak luminosity of high-redshift SNe la, therefore, these 
analyses have taken advantage of the empirical relation 
existing between the peak brightness and the light curve 
shape (LCS). Since this relation has been obtained from 
nearby SNe la only (Phillips 1993; Hamuy et al. 1995; 
Riess, Press & Kirshner 1995), it is important to examine 
whether it depends systematically on environmental prop- 
erties such as metallicity and age of the progenitor system. 
This Letter addresses the issue of whether a difference in 
the environmental properties is at the basis of the observed 
range of peak brightness. 



There are some observational indications that SNe la are 
affected by their environment. The most luminous SNe la 
seem to occur only in spiral galaxies, while both spiral and 
elliptical galaxies are hosts for dimmer SNe la. Thus the 
mean peak brightness is dimmer in ellipticals than in spi- 
ral galaxies (Hamuy et al. 1996). The SNe la rate per 
unit luminosity at the present epoch is almost twice as 
high in spirals as in ellipticals (Cappellaro et al. 1997). 
Moreover, Wang, Hoflich & Wheeler (1997) and Riess et 
al. (1999) found that the variation of the peak brightness 
for SNe located in the outer regions in galaxies is smaller. 

Hoflich, Wheeler, & Thielemann (1998) examined how 
the initial composition of the WD (metallicity and the 
C/0 ratio) affects the observed properties of SNe la. 
Umeda et al. (1999) obtained the C/0 ratio as a func- 
tion of the main-sequence mass and metallicity of the WD 
progenitors. In this Letter we suggest that the variation of 
the C/0 ratio is the main cause of the variation of SNe la 
brightness, with larger C/0 ratio yielding brighter SNe la 
(§ 2). We then show that the C/0 ratio depends indeed on 
environmental properties, such as the metallicity and age 
of the companion of the WD (§ 3), and that our model can 
explain most of the observational trends discussed above 
(§4). We then make some predictions about the brightness 
of SN la at higher redshift (§5). 
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2. EXPLOSION MODEL AND THE C/O RATIO OF WD 
PROGENITORS 

For the progenitors of SNe la, we adopt the single de- 
generate (SD) Chandrasekhar mass model, in which an 
accreting C-0 WD explodes when its mass reaches the 
critical mass Mja — 1.37— 1.38M0 (Nomoto, Thielemann, 
& Yokoi 1984). Merging of white dwarfs is likely to lead to 
accretion-induced-collapse rather than thermonuclear ex- 
plosion (Saio & Nomoto 1998). Chandrasekhar mass mod- 
els can reproduce well the spectrum and the light curves of 
SNe la, assuming either the explosion is induced by a de- 
flagration or by a delayed detonation (Hoflich & Khokhlov 
1996; Nugent et al. 1997). In these models, the brightness 
of SNe la is determined mainly by the mass of ^^Ni syn- 
thesized (MNise). Observational data suggest that MNise 
for most SNe la lies in the range Mnisg ~ 0.4 — 0.8Mq 
(e.g. Mazzali et al. 1998). This range of Mwise can result 
from differences in the C/O ratio in the progenitor WD as 
follows. 

In the deflagration model, a faster propagation of the 
convective deflagration wave results in a larger MNise- For 
example, a variation of the propagation speed by 15% in 
the W6 - W8 models results in MNise values ranging be- 
tween 0.5 and O.7Af0 (Nomoto et al. 1984), which could 
explain the observations. The actual propagation of the 
deflagration depends on the highly non-linear behavior of 
the turbulent flame (Niemeyer & Hillebrand 1995), and so 
it may be very sensitive to the C/O ratio. Qualitatively, a 
larger C/O ratio leads to the production of more nuclear 
energy and buoyancy force, thus leading to a faster propa- 
gation and a larger Mm56- Quantitatively, further studies 
of the turbulent flame are necessary to confirm that the 
expected range of C/O results in the required 15-20% vari- 
ation of the flame speed. 

In the delayed detonation model, MNise is predom- 
inantly determined by the deflagration-to-detonation- 
transition (DDT) density pddt, at which the initially 
subsonic deflagration turns into a supersonic detonation 
(Khokhlov 1991). We reproduce the relation between 
Pddt and MNise in Figure 1 by performing hydrodynam- 
ical calculations for several values of pddt, as in Nomoto 
et al. (1997), Kishimoto et al. (1999), and Iwamoto et 
al. (1999). The pre-explosive WD model is model C6 
(Nomoto et al. 1984), and the flame speed of the initial 
deflagration is assumed to be 3% of the local sound veloc- 
ity. Figure 1 shows that if the transition density varies in 
the range pddt — 1.3 — 3 x 10^ g cm~^, the resulting vari- 
ation of MNi56 is large enough to explain the observations. 

Possible mechanisms for DDT to occur have been stud- 
ied by Arnett & Livne (1994), Niemeyer & Woosley (1997), 
and Khokhlov, Oran & Wheeler (1997): when the defla- 
gration wave reaches a sufficiently low density, p ^ 10^ g 
cm~^, the turbulent motion associated with the flame may 
destroy the burning front. The resulting turbulent mix- 
ing between ashes and fuels efficiently heats up the fuel 
and could produce a region with a very shallow temper- 
ature gradient. In such a region, successive spontaneous 
ignitions cause the over-driven deflagration to propagate 
supersonically. This may induce a detonation wave if the 
mass of the region exceeds a critical mass AMddt- This 
critical mass is quite sensitive to the carbon mass fraction 
X(C), e.g. AMddt -^ 10"^^ and IQ-'^^Mq at p = 3 x 10^ 



g cm^'' for X{C) = 1.0 and 0.5, respectively (Niemeyer & 
Woosley 1997). Though the exact value of pddt is still de- 
bated, and its dependence on ^(C) has not been studied, 
it is not unlikely that for a larger X{C) DDT can occur at 
larger pddt. Hence a larger X{C) is likely to result in a 
larger pddt and Mnisg- 

In this Letter, therefore, we postulate that Mnisg and 
consequently brightness of a SN la increase as the progen- 
itors' C/O ratio increases, as illustrated in Figure 1, where 
the range of MNise ~ 0.5 — O.SMq is the result of an X{C) 
range 0.35 — 0.5, which is the range of ^(C) values of our 
progenitor models described below (Figure 2). The X(C) 
~ AfNi56 relation we adopt is still only a working hypothe- 
sis, which needs to be proved from studies of the turbulent 
flame during explosion. 

Hoflich et al. (1998) considered the dependence on 
X{C) but they assumed in their DDT model that pddt 
and ^(C) are independent parameters. They showed that 
for the same pddt a smaller X{C) leads to a slightly 
brighter SN la despite a slightly smaller Mnisg produced, 
because a smaller fraction of the explosion energy goes into 
the kinetic energy. In this Letter we assume that ^(C) is 
the primary parameter to determine pddt and thus Mjsh^q. 
The assumed variation of Mnisg is so large that a smaller 
X{C) yields an intrinsically dimmer SN. 



3. METALLICITY AND AGE EFFECTS 

In this section we discuss how the C/O ratio in the WD 
depends on the metallicity and age of the binary system. 
The C/O ratio in C+0 WDs depends primarily on the 
main-sequence mass of the WD progenitor and on metal- 
licity. According to the evolutionary calculations for 3—9 
Mq stars by Umeda et al. (1999), the C/O ratio and its 
distribution are determined in the following evolutionary 
stages of the close binary. 

1) At the end of central He burning in the 3—9 Mq pri- 
mary star, C/0< 1 in the convective core. The mass of 
the core is larger for more massive stars. 2) After central 
He exhaustion, the outer C-l-0 layer grows via He shell 
burning, where C/0> 1 (Umeda et al. 1999). 3a) If 
the primary star becomes a red giant (case C evolution; 
e.g. van den Heuvel 1994), it then undergoes the second 
dredge-up, forming a thin He layer, and enters the AGB 
phase. The C-l-0 core mass, Mco, at this phase is larger 
for more massive stars. For a larger Mco the total car- 
bon mass fraction is smaller. 4a) When it enters the AGB 
phase, the star greatly expands and is assumed here to 
undergo Roche lobe overflow (or a super-wind phase) and 
to form a C-l-0 WD. Thus the initial mass of the WD, 

Myyjj, in the close binary at the beginning of mass accre- 
tion is approximately equal to Mco- 3b) If the primary 
star becomes a He star (case BB evolution), the second 
dredge-up in (3a) corresponds to the expansion of the He 
envelope. 4b) The ensuing Roche lobe overflow again leads 

to a white dwarf of mass M^j-, = Mco ■ 

5) After the onset of mass accretion, the WD mass grows 
through steady H burning and weak He shell flashes, as de- 
scribed in the WD wind model (Hachisu, Kato, & Nomoto 
1996, 1999, and Hachisu et al. 1999; hereafter, HKN96, 
HKN99, and HKNU99, respectively). The composition 
of the growing C-l-0 layer is assumed to be C/0=1. 6) 



The WD grows in mass and ignites carbon when its mass 
reaches Mja — 1.367Mq, as in the model C6 of Nomoto 
et ah (1984). Because of strong electron-degeneracy, car- 
bon burning is unstable and grows into a deflagration for 
a central temperature of 8 x 10^ K and a central density of 
1.47 X 10^ g cm~^. At this stage, the convective core ex- 
tends to Mr = 1.14M0 and the material is mixed almost 
uniformly, as in the C6 model. 

In Figure 2 we show the carbon mass fraction X{C) in 
the convective core of this pre-explosive WD, as a function 
of metallicity (Z) and initial mass of the WD before the 
onset of mass accretion, Mco- Figure 2 reveals that: 1) 
X{C) is smaller for larger Mqo- 2) The dependence of 
X{C) on metallicity is small when plotted against Mco, 
even though the relation between Mco and the initial stel- 
lar mass depends sensitively on Z (Umeda et al. 1999). 

Metallicity dependent wind during mass accretion: In 
the SD Chandrasekhar mass model for SNe la, a WD ex- 
plodes as a SN la only when its rate of the mass accretion 
(M) is in a certain narrow range (e.g., Nomoto & Kondo 
1991). IIKN96 showed that the accreting WD blows a 
strong wind if M exceeds the rate Mb at which steady 
burning can process the accreted hydrogen into He. If the 
wind is sufficiently strong (i.e., the wind velocity v^ ex- 
ceeds the escape velocity Wosc of the WD), the WD can 
avoid the formation of a common envelope and increase 
its mass continuously at a rate Mb by blowing the extra 
mass away in a wind. 

In this model, which is adopted in the present study, an 
interesting metallicity effect has been found (Kobayashi 
et al. 1998; Hachisu & Kato 1999). The wind velocity is 
higher for larger Mwd and larger Fe/H because of higher 
luminosity and larger opacity, respectively. In order to 
blow sufficiently strong wind (i.e., v^ > Wcsc), Mwd should 
exceed a certain mass M^ (Fig. 6 of HKN99). As seen in 
Figure 1 of Kobayashi et al. (1998), M^ is larger for 
lower metalhcity; e.g., M^ = 0.65, 0.85, and 0.95 Mq for 
Z = 0.02, 0.01, and 0.004, respectively. In order for a 
WD to grow its mass at M > Mb, its initial mass Mco 
should exceed M^. In other words, M^ is the metallicity- 
dependent minimum Mco required for a WD to become an 
SN la {strong wind condition in Fig. 2). The upper bound 
Mco — 1.O7M0 is imposed by the condition that carbon 
should not ignite and is almost independent of metallicity. 

As shown in Figure 2, the range of Mco can be con- 
verted into a range of ^(C). From this we find the fol- 
lowing metallicity dependence for X(C): 1) The upper 
bound of X(C), which is determined by the lower limit 
on Mco imposed by the metallicity-dependent conditions 
for a strong wind, e.g., ^(C) < 0.51, 0.46 and 0.41, for 
Z=0.02, 0.01, and 0.004, respectively. 2) On the other 
hand, the lower bound, X{C) ~ 0.35 — 0.33, does not de- 
pend much on Z, since it is imposed by the maximum 
Mco- 3) Assuming the relation between Mwise and X{C) 
given in Figures 1 and 2, our model predicts the absence 
of brighter SNe la in lower metallicity environment. 

Age effects: In our model, the age of the progenitor sys- 
tem also constrains the range of ^(C) in SNe la. In the 
SD scenario, the lifetime of the binary system is essen- 
tially the main-sequence lifetime of the companion star, 
which depends on its initial mass M2. IIKNU99 and 
HKN99 have obtained a constraint on M2 by calculating 



the evolution of accreting WDs for a set of initial masses 
of the WD (M^^ ~ Mco) and of the companion (M2), 
and the initial binary period (Po)- In order for the WD 
mass to reach Mia, the donor star should transfer enough 
material at the appropriate accretion rates. The donors 
of successful cases are divided into two categories: one 
is composed of slightly evolved main-sequence stars with 
M2 ~ 1.7 - 3.6Mq (for Z=Om), and the other of red- 
giant stars with M2 ~ 0.8-3.1Mq (for Z=0.02) (HKN99, 
HKNU99; also Li & van den Heuvel 1997). 

If the progenitor system is older than 2 Gyr, it should be 
a system with a donor star of M2 < I.7M0 in the red-giant 
branch. Systems with M2 > I.7M0 become SNe la in a 
time shorter than 2 Gyr. Likewise, for a given age of the 
system, M2 must be smaller than a limiting mass. This 
constraint on M2 can be translated into the presence of a 
minimum Mco for a given age, as follows: For a smaller 
M2, i.e. for the older system, the total mass which can 
be transferred from the donor to the WD is smaller. In 
order for Mwd to reach Mja, therefore, the initial mass 
of the WD, M^\^ ~ Mco, should be larger. This im- 
plies that the older system should have larger minimum 
Mco as indicated in Figure 2. Using the Ar(C)— Mco and 
MNi56— X(C) relations (Figs. 1 and 2), we conclude that 
WDs in older progenitor systems have a smaller X[C), and 
thus produce dimmer SNe la. 

4. COMPARISON WITH OBSERVATIONS 

The first observational indication which can be com- 
pared with our model is the possible dependence of the SN 
brightness on the morphology of the host galaxies. Hamuy 
et al. (1996) found that the most luminous SNe la occur in 
spiral galaxies, while both spiral and elliptical galaxies are 
hosts to dimmer SNe la. Hence, the mean peak brightness 
is lower in elliptical than in spiral galaxies. 

In our model, this property is simply understood as the 
effect of the different age of the companion. In spiral galax- 
ies, star formation occurs continuously up to the present 
time. Hence, both WD+MS and WD+RG systems can 
produce SNe la. In elliptical galaxies, on the other hand, 
star formation has long ended, typically more than 10 Gyr 
ago. Hence, WD-I-MS systems can no longer produce SNe 
la. In Figure 3 we show the frequency of the expected SN I 
for a galaxy of mass 2 x IO^Mq for WD-hMS and WD-FRG 
systems separately as a function of Mco ■ Here we use the 
resuhs of HKN99 and HKNU99, and the Mco - ^(C) 
and MNi56 — -'^(C) relations given in Figure 2. Since a 
WD with smaller Mco is assumed to produce a brighter 
SN la (larger Mnisg), our model predicts that dimmer SNe 
la occur both in spirals and in ellipticals, while brighter 
ones occur only in spirals. The mean brightness is smaller 
for ellipticals and the total SN la rate per unit luminosity 
is larger in spirals than in ellipticals. These properties are 
consistent with observations. 

The second observational suggestion is the radial distri- 
bution of SNe la in galaxies. Wang et al. (1997) and Riess 
et al. (1998) found that the variation of the peak bright- 
ness for SNe la located in the outer regions in galaxies is 
smaller. This behavior can be understood as the effect of 
metallicity. As shown in Figure 2, even when the progen- 
itor age is the same, the minimum Mco is larger for a 
smaller metallicity because of the metallicity dependence 
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of the WD winds. Therefore, our model predicts that the 
maximum brightness of SNe la decreases as metallicity de- 
creases. Since the outer regions of galaxies are thought to 
have lower metallicities than the inner regions (Zaritsky, 
Kennicutt, Huchra 1994; Kobayashi & Arimoto 1999), our 
model is consistent with observations. Wang et al. (1997) 
also claimed that SNe la may be deficient in the bulges of 
spiral galaxies. This can be explained by the age effect, 
because the bulge consists of old population stars. 

5. CONCLUSIONS AND DISCUSSION 

We have suggested that ^(C) is the quantity very likely 
to cause the diversity in AfNise and thus in the brightness 
of SNe la. We have then shown that our model predicts 
that the brightness of SNe la depends on the environment, 
in a way which is qualitatively consistent with the observa- 
tions. Further studies of the propagation of the turbulent 
flame and the DDT are necessary in order to actually prove 
that ^(C) is the key parameter. 

Our model predicts that when the progenitors belong 
to an old population, or to a low metal environment, the 
number of very bright SNe la is small, so that the variation 
in brightness is also smaller. In spiral galaxies, the metal- 



licity is significantly smaller at redshifts z > 1, and thus 
both the mean brightness of SNe la and its range tend to 
be smaller. At z > 2 SNe la would not occur in spirals at 
all because the metallicity is too low. In elliptical galaxies, 
on the other hand, the metallicity at redshifts z ~ 1 — 3 
is not very different from the present value. However, the 
age of the galaxies at z ~ 1 is only about 5 Gyr, so that 
the mean brightness of SNe la and its range tend to be 
larger at z > 1 than in the present ellipticals because of 
the age effect. 

We note that the variation of X{C) is larger in metal- 
rich nearby spirals than in high redshift galaxies. There- 
fore, if X{C) is the main parameter responsible for the 
diversity of SNe la, and if the LCS method is confirmed 
by the nearby SNe la data, the LCS method can also be 
used to determine the absolute magnitude of high redshift 
SNe la. 
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FiG. 1. — Deflagration to detonation transition density (/9ddt) vs. ^^Ni mass (AfNise)- The ^(C) axis is approximately 
jiven as discussed in § 2. These results are obtained for X{C)=0A3. For a given pddt varying X{C) between X(C)=0.35 
- 0.50 changes the ^"^Ni mass only by -- ±0.02Mq. 
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Fig. 2. — The total ^^C mass fraction included in the convective core of mass, M — 1.14M0, just before the SN la 
explosion as a function of the C+0 core mass before the onset of mass accretion, Mco- The lower bounds of Mqo 
obtained from the age effects and the conditions for strong wind to blow are also shown by arrows. The axis of M^isg is 
obtained from the X{C) — MNise relation assumed in Figure 1. 
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Fig. 3. — SN la frequency for a galaxy of mass 2 x 10^^ Mq as a function of Mco for Z=0.02. For the WD+RG system, 
constraints from the companion's age are shown by the arrows. SNe la from the WD+MS system occur in spirals but 
not in ellipticals because of the age effect. Mco and MNise can be related as shown here if the X{C) — Mnisg relation in 
Figure 1 is adopted. 



